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Considering the supersymmetric models with light neutralino and A-parity violation, we perform 
estimates of the signal rate expected at the recently proposed fixed target SHiP experiment exploiting 
the CERN SPS beam of 400 GeV protons. We extend the existing studies by introducing new 
production channels (in particular through the beauty mesons) and decay modes. We also constrain 
the model parameter space from analysis of negative results of the CHARM experiment. 


I. INTRODUCTION 

Supersymmetric (SUSY) extensions of the standard 
model of particle physics (SM) provide taming of the ra¬ 
diative corrections to the Higgs boson mass (for a review 
see [1]). The crucial role is played by the superpartners of 
the SM particles, whose contributions cancel the quadrat- 
ically divergent quantum corrections of the SM particles. 
Therefore, one naturally expects to observe (some of) the 
superpartners at the TeV mass scale. Testing this pre¬ 
diction is one of the main tasks of run 2 at the LHC. 

However, there can be renegades in the SUSY world 
with masses (much) below the TeV scale and suffi¬ 
ciently suppressed interactions with the SM fields, so they 
are missed by numerous previous searches and are be¬ 
yond the reach of the LHC. The renegade hunting can 
be performed at a beam-dump experiment, where the 
high statistics of proton(electron)-proton collisions might 
compensate for smallness of the couplings, so the new 
light particles can be produced. Recently a proposal has 
been submitted [2, 3] to build the new experiment at 
CERN with a 400 GeV SPS proton beam. Originally mo¬ 
tivated as a facility to search for sterile neutrinos (heavy 
neutral leptons) of 0(1) GeV mass [2, 4, 5], lately it has 
been recognized as a universal tool to probe various mod¬ 
els predicting light, sufficiently long-lived, neutral parti¬ 
cles; it has been named SHiP (Search for Hidden Parti¬ 
cles) [3]. The SHiP physics case is presented in a separate 
paper [6], including a number of the SUSY renegades. 

In this work, we consider the light unstable neutralino 
in supersymmetric models with R-parity violation (for 
reviews, see [7 9]). R-parity is a discrete multiplicative 
symmetry ascribing factor, 

R p = (~1) 3B+L+2S , (1) 

to any particle of baryon charge B , lepton charge L , and 
spin S. All SM fields (including scalars of the extended 
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Higgs sector) have R p = +1, while their superpartners 
have R p = —1. If R-parity is conserved, superpartners 
may only be created in pairs. R-parity guarantees stabil¬ 
ity of the lightest supersymmetric partner (LSP) which 
becomes a candidate to form the dark matter component 
of the Universe. 

However, the theoretical grounds of R-parity have been 
questioned (see, e.g., Refs [10, 11]), and one can introduce 
R-parity violating (RPV) terms in the Minimal Super- 
symmetric extension of the SM (MSSM). These terms 
trigger decays of superpartners into SM particles, in par¬ 
ticular, decays of the LSP. The latter may be the light¬ 
est mass eigenstate in the sector of neutral fermion su¬ 
perpartners called neutralinos. If the neutralino is suf¬ 
ficiently light, it is an example of the renegades to be 
searched for at the SHiP experiment. There, the neu¬ 
tralinos can be produced either directly in proton scat¬ 
tering off the target material or indirectly via decays of 
secondary hadrons, and they later decay into SM parti¬ 
cles exhibiting signatures very similar to those of sterile 
neutrinos [4]. Thus, the procedures applied in data anal¬ 
ysis to probe both models are very similar. The differ¬ 
ence between the two models is expected in the produc¬ 
tion channel pattern and in decay channel patterns, so 
that the momentum spectra of each final state and the 
weights of the different final states (A^/pF, etc) 

are generically not the same in the two models. 

Light neutralino phenomenology has previously been 
studied in the SHiP physics paper [6]. Here, we consid¬ 
erably extend this study by including more production 
(from neutral charm mesons and from beauty mesons) 
and decay (tt^P, V^V^v) channels and by obtaining lim¬ 
its on the model parameter space from analysis of the 
published results of the CHARM experiment [12, 13]. 

The paper is organized as follows. In Sec. II we intro¬ 
duce the model and calculate the neutralino production 
rates and decay rates for a number of channels. In Sec. Ill 
we give the estimate of the number of signal events ex¬ 
pected in the SHiP fiducial volume. In Sec. IV we present 
the sensitivity of the SHiP experiment to the model pa¬ 
rameters (assuming zero background) and find new limits 
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on the model parameters from the published results of the 
CHARM experiment. We summarize in Sec. V. 


II. SUPERSYMMETRY WITH RPV: 

PRODUCTION AND DECAYS OF LIGHT 
NEUTRALINOS 

R-parity is explicitly violated by the following terms 
in the MSSM superpotential 

W n „ = \ ijk e ah L1L)E^ + A u fc e a6 L?Q$I>£ (2) 

+X^ jk e aM Uf' a Df‘ 8 D^ + me ab L“H% , (3) 

where dimensionless couplings A ij k and mass parameters 
Pi ( i,j,k run over the three matter generations) charac¬ 
terize violation of R-parity, the superscript C refers to the 
charge conjugated fields, indices a, b = 1,2 indicate the 
SU(2)w doublet components (L and Q are lepton and 
quark doublets; E 1 D, and U are lepton , down-type, 
and up-type quark singlets, respectively), while a,/3,j 
count SU(3)c triplet components; e ab and e aj a 7 are fully 
antisymmetric 2x2 and 3x3x3 tensors. Now, if all the 
terms (2), (3) are present, they initiate the fast proton 
decay. This process can be forbidden with some discrete 
remnant of R-parity. In particular, the baryon triality 
[14] forbids the first term in (3) and hence keeps the 
proton stable. In what follows we concentrate on the 
phenomenology of the RPV terms in Eq. (2) and neglect 
the terms in Eq. (3). Whereas the lightest neutralino 
should be heavier than 46 GeV in the constrained MSSM 
with five parameters [15], the authors of Ref. [16] show 
that this bound can be relaxed and, even a massless neu¬ 
tralino is possible. In this study we consider models with 
the mass of the lightest neutralino in a GeV range. 


Neutralinos as the LSPs could be created in decays 
of heavier sparticles. Missing momentum carried out by 
the LSP, which escaped from a detector, remains one of 
the main signatures in collider searches of supersymme¬ 
try. The proposed center-of-mass energy of the SHiP 
experiment y/s = 27.4 GeV is too low to create heavy 
superpartners with masses above the electroweak scale, 
as we anticipate from LEP-II, Tevatron, and LHC run 
I. One can check that the neutralino direct production 
in the proton-proton scattering is negligibly low. There¬ 
fore, here we study an indirect production of neutralinos 
in decays of heavy mesons via R-odd couplings A' in (2). 
These R-odd couplings also lead to neutralino decays into 
the ordinary SM particles. 

A. Neutralino production in decays of heavy 
mesons 


Light enough neutralinos can be produced in decays 
of heavy mesons (charm D and beauty B) provided that 
we have R-parity-violating coupling A 7 as shown in Fig. 1. 



FIG. 1. Typical Feynman graphs of neutralino production in 
meson decay via A-parity-violating couplings A 7 . 

Expressions for the partial widths of B 0 and B + meson 
decays can be found in Ref. [17] (for details and derivation 
of similar expressions see also Ref. [18]). For a neutralino 
being a pure bino state, they read 
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where p cm is the 3-momentum of outgoing particles in 
the rest frame of decaying meson, m u , md , m b are quark 
masses, me t is mass of the final state lepton (electron 
or muon), Mp, , M^ r , ...are sfermion masses, Y Vi , 
Yd, Y b , ... are corresponding hypercharges (coming from 
neutralino couplings in the pure bino limit), g' is U(l)y 
gauge coupling constant, M^o is neutralino mass, Mg+, 
M b o are B + and B° masses, respectively, and f B = 204± 


30 MeV [15] is the R-meson decay constant. Note that 
definition of the constant f B in [17] differs from that by 
the Particle Data Group [15]. Therefore, formulas from 
[17] should be multiplied by 1/2 for our choice of f B . 

Without loss of generality, hereafter we assume the 
common mass scale of sfermions Mj = Mp t = = 

■ ■ ■ = Mj jR This assumption simplifies further phe¬ 
nomenological treatment since Eqs. (4) and (5) turn into: 
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and the production rates of neutralinos are proportional 
to the squared combination A' 13 /M?. Note that our as¬ 
sumption is not related to any particular pattern of mass 
spectrum of RPV SUSY. In a more general case of differ¬ 
ent masses of superpartners, the dominant contribution 
will come from the intermediate sfermion with the high¬ 
est value of X'/Mj, e.g., from the lightest sfermion, if all 
RPV couplings are of the same order. 

Expressions similar to (6) and (7) can be obtained for 
decays of D mesons with obvious replacement of cou¬ 
plings A' 13 —> A' 21 and replacement of the corresponding 
quark flavors. 


B. Decay pattern 


three-body leptonic decay via Xijk, two-body semilep- 
tonic decay into a charged pion via the coupling A' ll7 
and two-body semileptonic decay into a charged kaon via 
A' 12 . These states yield the detectable at SHiP signature 
of two charged particles originated from a single vertex. 

The amplitude of the leptonic neutralino decay 
through a virtual slepton or sneutrino was calculated in 
Ref. [19]. The decay width in the pure bino limit (ne¬ 
glecting masses of the final-state particles) reads 


r(x? 
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We study three different decay channels of light neu¬ 
tralinos with two charged particles in the final state: 


The rate of semileptonic decay can be calculated analo¬ 
gously to Eq. (4). In the pure bino limit one has [6] 


T(x? -> K-ef) 
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and the expression analogous to (9) for 
which is proportional to |A' n | 2 . One can see from 
Eqs. (8) and (9) that the decay rates of a neutralino, as 
well as its production rates [Eqs. (6) and (7)], are pro¬ 


portional to the factor 



If Xi is produced in decays of D mesons via A' 21 then 


there is one additional decay channel Xi ~t K °i/j. A 
special study is needed to understand whether it can be 
distinguished from background since there are no charged 
particles in the final state [20] but it affects the neutralino 
lifetime. Hence we take into account the decay rate 


T(X? K%) 
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( 10 ) 


and a similar one for Xi 7r °^i (which is proportional III. SIGNAL EVENT RATES 

to |A' 1:l | 2 ) in further investigations. 

The goal of the present paper is to estimate the event 
rate of light neutralino decays within the SHiP detec¬ 
tor. We consider the light neutralinos created in decays 
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of heavy mesons. These heavy mesons are produced, in 
turn, by 400 GeV protons scattering off the target mate¬ 
rial. 

For the production cross section of a particle (a neu- 
tralino in our case) with 3-momentum p created in a de¬ 
cay of the heavy meson H (H = D , B) with 3-momentunr 
k. one has: 


The probability of Xi decay inside the fiducial volume 
of the SHiP detector is 


— ish/Eo /i — ilid/bo \ ^fid 

Wdet = e *i(l - e x i) ~ —, 

1 p 

Xl MyoT 

Xi 


(15) 


= e /* ,(Si) iSk' <u) 

where p = \p\,k = |fc|, B is the branching ratio of H 
two-body decay to neutralino, f(p , k) is the momentum 
distribution of a neutralino, and dkd g// d( j }k the differen¬ 
tial production cross section of meson H in pp collisions. 
All the momenta in Eq. (11) are given in the labora¬ 
tory frame. In the rest frame of the meson H (denoted 
by an asterisk), neutralino 3-momentum is uniformly dis¬ 
tributed and one finds: 

f(p*,0) = ^6(p*»Pl-m 2 ), (12) 

where p * At is the 4-momentum of the neutralino in the 
rest frame of the decaying meson. Note that the value of 
p* = | p* | is fixed by kinematics of the two-body decay. 
In order to boost expression (12) to the laboratory frame, 
one should multiply it by the appropriate Jacobian and 
express p* = p* ( p ,, k ) in terms of the 3-monrenta of the 
neutralino (p) and the decaying meson ( k ) in the labora¬ 
tory frame. 

The differential cross section of D-meson production 
in pp interactions at the center-of-mass energy y/s = 
27.4 GeV, which is relevant for the SHiP setup, was mea¬ 
sured by the LEBS-EHC Collaboration [21]. It has been 
found that the differential production cross section is well 
represented by the empirical form [21] 

d ^p 2 T = \ ( a ( D /D)(n + 1)6) (1 - \x F \) n exp {-bp%) 

(13) 

with n = 4.9 ± 0.5, b = (1.0 ± 0.1) GeV . We adopt the 
same value of the inclusive D/D cross section <j(D/D ) = 
18pb [6]. The differential cross section in (13) depends 
on transverse pr and longitudinal pl components of 3- 
momenta through x F = 2 Pl/\[s and can be related to 
that used in Eq. (11) as follows: 

d 3 UD Ak 2 sm0 k d 3 a 

dkd9 k d(j) k ^/~s dxpdk^d(f>k 

The beauty production cross section has not been mea¬ 
sured in the interesting energy region. Therefore, fol¬ 
lowing the SHiP Collaboration [6] we extrapolate ex¬ 
isting data [22] to estimate the number of produced B 
mesons. To estimate the angular distribution of pro¬ 
duced B mesons, we employ the theoretical results from 
Ref. [23] together with the Lund fragmentation model 
[24]. 


with Z s h denoting the muon shielding length (the distance 
between the collision point and the detector, 63.8 m for 
SHiP [3]) and lad referring to the length of the detector 
fiducial volume (60 m); the second equality in (15) is valid 
when Is h < J*o. 

The proposed geometry of the SHiP detector [3] is a 
60 nr length cylindrical vacuum tank with an elliptical 
section of x and y semiaxes 2.5 m and 5 nr long, respec¬ 
tively. In the following estimates of the signal event num¬ 
bers, we utilize a more conservative fiducial volume that 
is a cone formed by the vertex in the target and the 
5 nr x 10 m ellipse at the very end of the fiducial volume. 
It covers part of the elliptical section. We argue that this 
choice is quite reasonable since neutralino decay prod¬ 
ucts should be tracked by the detector placed at the end 
of the vacuum tank. We select only neutralinos x r with 
3-nromenta inside the cone region described above. 

The number of neutralino decays within the detector 
is given by 


N = 


iVpoT 

^pp,total 



Wdet 


da^o 

_ A-l 

dpdOdcj) 


d 3 P 


(16) 


where the distribution of neutralinos over the 3- 
momentunr is defined in (11) and “cut” refers to the con¬ 
straint on the neutralinos’ 3-momenta described above 
and IVpoT = 2 x 10 20 is the number of protons on target 
during 5 years of operation [6]. 

We consider the four production channels of light neu¬ 
tralinos described in Sec. II: D ± , D° decays via coupling 
X', m and B ± , B° decays via coupling A' 31 . According 
to the SHiP technical proposal [3], at least two charged 
particles are required to distinguish the signal from a 
background. Hence we consider three decay channels: 
Xi -t ttvjtk, Xi Xi 7 t~1+, driven by 

\jki Xl2> •Mil correspondingly. As a result, we have 
six combinations of couplings that can be tested by SHiP. 

Resulting event rates for neutralinos produced in D° 
and B° decays are shown in Fig. 2. Plots for the D + and 
B + channels are very similar to those in Fig. 2. At small 
couplings the number of events drops due to poor pro¬ 
duction, while at relatively large couplings the number of 
events also drops because of fast neutralino decay; thus, 
one has L-o <C l s h and the neutralino flux in the detector 
is exponentially suppressed [see Eq. (15)]. The num¬ 
ber of events depends both on combination (A/Af~) 2 and 
on neutralino mass. To demonstrate this dependence we 
solve Eq. (16) with N = 3 and find ( X/Mj ) 2 as a func¬ 
tion of A/^o. As one can see from Fig. 2, there are two 
solutions to this equations: one corresponding to small 
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FIG. 2. Number of signal events for neutralinos produced in 
D° decay (upper panel) and in B° decay (lower panel). De¬ 
cays of the neutralino to e + Vie~ (solid line), to n + £~ (dashed 
line), and to K + l~ (dotted line) are considered. A stands for 
the appropriate RPV couplings. The horizontal line repre¬ 
sents three events: their absence implies an upper limit on 
the model parameters placed at 95% confidence level within 
the Poisson statistics. 


couplings and slow decay and another corresponding to 
relatively large couplings and fast decay. Since large cou¬ 
plings are excluded by different searches (see, e.g., Ref. 
[7] for details), we show in Fig. 3 a solution corresponding 
to small couplings. Note that for the small couplings the 
probability of neutralino decay inside the fiducial volume 
(15) and, consequently, the number of events is propor¬ 
tional to the decay width. Therefore, one can use Eqs. 
(4) - (7) and Eqs. (8) - (10) (see also expressions from 
Ref. [19]) in order to rescale X'/Mj dependence for the 
generic case of a nondegenerate mass spectrum and any 
pattern of RPV couplings. 

One can see from Fig. 3 that the decay channel x? —t 
K°i/j sufficiently affects the lifetime and, subsequently, 
the event rate of neutralinos created in D-meson decays 
via the coupling X' i21 . 


IV. SHIP SENSITIVITY TO AND CHARM 
BOUNDS ON RPV 






FIG. 3. Expected sensitivity of the SHiP experiment to a light 
neutralino with RPV. The solid line refers to leptonic decay, 
the dashed line is for decay to e + 7r“, and the dotted line is 
for decay to e + K~. A stands for the appropriate combination 
of RPV couplings. 


Absence of the events, while the three signal events are 
expected, implies 95% confidence level bounds (if back¬ 


ground is negligible, which is true in our case [3]). Lim¬ 
its on RPV couplings that could be placed by the SHiP 
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experiment depend both on the common sfermion mass 
scale and on the neutralino mass. Exclusion limits on var¬ 
ious combinations of RPV couplings are shown in Fig. 4. 
Dark shaded regions are excluded by previous stud¬ 
ies. Namely, the constraints on A' 121 , A' 113 , Am, A 122 , A 123 
have been obtained from charged current universality, 
and the constraint on A 311 has been obtained from neutri¬ 
noless double beta decay (see Refs. [25], [7] for details). 
Note that these constraints scale as oc X/Mj, while in our 
case of the SHiP sensitivity the scaling goes as oc A/M?. 

Vertical lines in Fig. 4 at Mr = 1 TeV show the mass 
scale that is excluded by the LHC searches in the case 
of approximately equal sfermion and gluino masses [26]. 
Note, however, that for some particular RPV spectra this 
bound can be as low as 800 GeV [27]. 

We list our estimates of SHiP bounds on these cou¬ 
plings in Table I. As one can see from Fig. 3 the bounds 

TABLE I. Estimates of SHiP sensitivity to and CHARM 

bounds on combinations of RPV couplings. In the first three 

rows we set AUo = 1 GeV and AUo = 4 GeV for the last three 
Xi x 1 

rows. Indices j, k = 1,2 and i = 1, 2, 3 indicate flavor of the 
final-state leptons. 


A 

Expected sensitivity 
SHiP, Mf/ TeV 1 2 3 

Upper limit 
CHARM, AL 2 /TeV 2 

\J Al 21 A ijk 

2.4 x 10 -3 

2.5 x 10~ 2 

\J X'\2lX'jU 

1.2 x 10 -3 

- 

yj A121A/2I 

1.4 x 10 -3 

- 

co 

> 

2.4 x 10 -3 

2.5 x 10~ 2 

yA' 113 A ' n 

3.9 x 10 -3 

- 

^/a' 113 a ' 21 

4.0 x 10 -3 

- 


listed in Table I are valid for a wide range of kinemati¬ 
cally allowed region (with phase-space corrections at the 
boundaries) of M^o. 

In order to illustrate advantages of the SHiP facility, 
we also present in Fig. 4 our estimates of the bounds on 
RPV couplings that follow from the absence of signal in 
the CHARM experiment. These bounds are obtained by 


adapting the whole procedure described above for the 
CHARM geometry [12]. The CHARM experiment has 
exploited the same 400 GeV beam as the SHiP plans. 
The detector was located at 480 m downstream from the 
beam dump. Therefore, it covers a sufficiently smaller 
solid angle compared to the SHiP. The length of the de¬ 
cay region was 35 m and the radius of the calorimeter 
placed at the end of decay volume was 1.5 m. The total 
amount of protons on target equalled 2.4 x 10 18 [12]. To 
the best of our knowledge there was no special investiga¬ 
tion of the CHARM sensitivity to RPV SUSY, but very 
similar signatures of heavy neutral lepton decays to the 
SM leptons were studied in Refs. [12], [13]. 

V. CONCLUSION 

To summarize, we have estimated the sensitivity of the 
recently proposed SHiP experiment to the supersymmet¬ 
ric extensions of the SM with light neutralinos and 77- 
parity violation. For the i?-parity violating couplings A 
of order one, the SHiP will allow us to probe the super¬ 
partner mass scale as high as 30 TeV (see Fig. 3), which 
is in agreement with previous estimates in Ref. [ 6 ]. The 
number of signal events scales as oc (A'/M?) 4 . As a by¬ 
product we have obtained limits on the model parameters 
from nonobservation of anomalous events in the CHARM 
experiment. With respect to the CHARM, the SHiP will 
improve the sensitivity to R-parity-violating couplings by 
an order of magnitude. 

Several remarks are in order. First, other final states 
such as the mentioned neutral kaons must be considered 
as well. Second, the light neutralinos can be produced in 
pairs due to R-parity-conserving couplings, and the cor¬ 
responding new production channels can also be studied. 
Third, not only heavy meson but also heavy baryons can 
decay into light neutralinos, which gives additional pro¬ 
duction channels. Fourth, secondary hadrons, produced 
in the hadron showers initiated by 400 GeV proton scat¬ 
tering off target materials, can contribute to the light 
neutralino production. Fifth, r leptons produced mostly 
in decays of D s mesons allow us to probe other R-parity- 
violating couplings A,^, which is also worth investigat¬ 
ing. 

This work has been supported by Russian Science 
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FIG. 4. Bounds on A as a function of sfermion mass for the 
light neutralino mass M^o = 1 GeV. The region above the 
thick solid black line has been excluded by previous studies 
[25], [7]. The region above the gray line could be excluded by 
SHiP data. The dashed gray line corresponds to our estimate 
of CHARM bounds. The region to the left of the vertical line 
at Mj = 1 TeV is generally disfavored as the superpartner 
mass scale due to searches at the LHC (see discussion in the 
main text). 







